Retardation effects and the Coulomb pseudopotential in the theory of 

superconductivity 
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In the theory of electron-phonon superconductivity both the magnitude of the electron-phonon 
coupling A as well as the Coulomb pseudopotential /i* are important to determine the transition tem- 
perature Tc and other properties. We calculate corrections to the conventional result for the Coulomb 
pseudopotential. Our calculation are based on the Hubbard-Holstein model, where electron-electron 
and electron-phonon interactions are local. We develop a perturbation expansion, which accounts for 
the important renormalization effects for the electrons, the phonons, and the electron-phonon ver- 
tex. We show that retardation effects are still operative for higher order corrections, but less efficient 
due to a reduction of the effective bandwidth. This can lead to larger values of the pseudopotential 
and reduced values of Tc- The conclusions from the perturbative calculations are corroborated up 
to intermediate couplings by comparison with non-perturbative dynamical mean-field results. 

PACS numbers: 71.10.Fd, 71.27.+a,71.30.+h,75.20.-g, 71.10.Ay 



I. INTRODUCTION 

More than a century after its discovery superconduc- 
tivity continues to be subject of intense research in con- 
densed matter physics. The research ranges from topics 
geared to the technical application of the phenomenon 
over the search for new superconducting materials to fun- 
damental questions of microscopic mechanisms fi The lat- 
ter are important as one might hope that a good under- 
standing of the mechanisms will make the design of new 
superconductors at elevated temperatures feasible. ^'^ 
There are numerous classes of materials for which differ- 
ent mechanisms are discussed.— In most cases a bosonic 
pairing function is invoked, which could be of purely elec- 
tronic origin such as spin fluctuations, or the conventional 
phonon mechanism, but also more exotic mechanisms 
have been proposed. The electron-phonon mechanism 
has the longest history and probably the most estab- 
lished mathematical foundation. An effective electron- 
electron attraction generated by the electron-phonon 
interaction is part of the celebrated Bardeen-Cooper- 
Schrieffer (BCS) theoryj^ The more elaborate theory in- 
cluding microscopic details goes under the name Migdal- 
Eliashberg (ME) theory^ With the help of ME theory the 
superconducting properties of many elements and numer- 
ous alloys have been described accurately.^'* 

The general ideas of ME theory can be presented in 
relatively simple fashion, although the details of the 
complete framework, its foundations and specific appli- 
cations, involve a substantial degree of sophistication 
to which many researcher have contributed over the 
yearsii^"-t2 One cornerstone is Migdal's theorem, ^-"^ which 
employs the fact that the typical electronic energy scale 
Eel and the typical phonon energy scale cjph differ largely, 
such that Ec\/ijJph is of the order 100 and larger - the elec- 
trons move much faster than the phonons. It can then 
be shown that the perturbation theory of the electron- 



phonon problem greatly simplifies since vertex correc- 
tions are small. The infiuence of the bosonic pairing 
function a'^F{uj) on the electronic properties and the oc- 
currence of superconductivity can therefore be computed 
reliably. A remarkable aspect of the ME theory is that 
it is even well justified for large values of the electron- 
phonon coupling parameter A > 1 as long as the effec- 
tive expansion parameter ^ Awph/i^ci remains smallfii"— 
With the help of the ME equations, the pairing func- 
tion a'^F{Lu) and the phenomenological parameter for the 
Coulomb pseudopotential ^* were extracted from tunnel- 
ing measurement in a procedure termed inverse tunnel- 
ing spectroscopy!^ Based on those other properties such 
as Tc or thermodynamic quantities can be computed. 
The agreement with the respective experimental values of 
many elements and alloys, notably Pb and Nb, is within 
the range of a few percent.^ This leads to a consistent pic- 
ture and was taken as proof of the validity of ME theory 
and the electron-phonon mechanism/'* Further support 
comes from first principle calculations for a^F^uj), which 
are in many cases in remarkable agreement with the re- 
sults extracted from tunnelingf^i^ d^d^ 

In addition to the effective attraction induced by the 
electron-phonon coupling the applications of ME theory 
include the effects of the Coulomb repulsion. Due to 
the enormous success of the theory it is sometimes un- 
derstated that in contrast to the electron-phonon prob- 
lem no rigorous arguments exist for the treatment of the 
Coulomb repulsion,^*' the reason being the absence of a 
small parameter as in Migdal's theorem. Traditionally, 
this is seen as minor deficiency based on the following ar- 
guments. Some effects of the Coulomb interaction, such 
as the renormalization of the electron-phonon couplings, 
g, the electron and phonon dispersion Efc, are implicit 
in the experimentally derived results or in the first princi- 
ple calculations. It remains to deal with the direct repul- 
sion in the pairing channel which opposes s-wave super- 
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conductivity. Morel and Anderson proposed a procedure 
in two stagesf^^ first the Coulomb interaction is screened 
and averaged over the Fermi surface. In a second step 
it is projected to the phonon scale. The result possesses 
the famous iorm.^^^ 



(1) 



which is sometimes termed the Morel-Anderson (MA) 
pseudopotential. Here, fic is a dimensionless quantity 
which consists of a product of the the averaged, screened 
Coulomb interaction with the density of states per spin 
at the Fermi energy; Ed is an electronic scale, such as 
the half bandwidth D or the Fermi energy Ep, and Wph 
the phonon scale, e.g. the Debye frequency. Eq. (|T]) 
has the important property that for the typical large 
energy separation between electronic and phonon scale. 
Eel ^ Wph, such that log(i?oi/wph) ^5—10, /i* is much 
smaller than /ic. This effect is remarkable, as it enables 
the electron-phonon s-wave superconductivity to be pos- 
sible in spite of the Coulomb repulsion, which on a bare 
level is much larger and working directly against it. This 
is sometimes even termed the true mechanism of electron- 
phonon superconductivity*^ The textbook physical pic- 
ture is that electrons do not need to be close in position 
space and suffer from the Coulomb repulsion in order to 
pair since the electron-phonon interaction is retarded and 
thus electrons can pair with a "time-delay". Using appro- 
priate energy scales in Eq. ([T]) leads to estimates of the 
order /i* ^ 0.1 — 0.14. This fits well to the parameters fj,* 
obtained from tunneling in many elemental superconduc- 
tors and alloySfi Hence, the Coulomb effects are generally 
considered to be well described by fj,* , which is a fairly 
universal quantity. In contrast to the pairing function 
it is usually not calculated based on first principles and 
rather used as a fitting parameter. It is worth mention- 
ing that there is also an alternative ab-initio approach to 
superconductivity based on the DFT framework. ^%"= 

Review of the literature however also reveals evidence 
that the description of the Coulomb repulsion in the form 
of Eq. ([1]) is incomplete. In some cases, such as V or 
NbsGe,^ the values for /i* appear to be of the order 
0.2 — 0.3 substantially larger than the traditional quotes, 
even though the maximal phonon scales does not seem to 
be particularly large. Density functional theory (DFT) 
calculatior^i^ for a^F^uj) find good agreement with the 
tunneling results, but to explain the experimental val- 
ues for Tc also somewhat larger values of /z* have to be 
used. A specific case which raises doubts about the con- 
ventional framework is the example of elemental lithium 
at ambient pressurCi ^^i^^ The coupling constant was es- 
timated to be A ~ Q 4 ^22^24 -^Yhich seems to be in line 
with specific heat measurements. With the usual value 
of /i* ^ 0.1 and the appropriate phonon scale this implies 
Tc ~ IK. This is in contrast to experiments where for a 
long time no superconductivity was observed down to val- 
ues of 6mK)^ and only recently, Tc ^ 0.4mK was found 
at ambient pressure, which requires /i* ~ 0.23.^^ In this 



context, we also mention alkali doped picene, which was 
recently discovered to be superconducting at Tc ^ 7K 
and Tc ^ 18K depending on preparation..^- First prin- 
ciple calculations only seem to be able to explain these 
values of Tc based on an electron-phonon mechanism with 
relatively large values of fi* ~ 0.23.'^^ However, also other 
interpretations exist i^i^S 

It is our endeavor to reanalyze the expression for the 
Coulomb pseudopotential in Eq. ([T]) and calculate cor- 
rections to it. We focus on the reduction of phonon 
induced s-wave superconductivity due to the Coulomb 
repulsion between electrons. Superconductivity which 
is induced in an anisotropic higher order angular mo- 
mentum channel by purely repulsive interactions, such 
as the well-known Kohn-Luttinger effecli^, is not dealt 
with in the present work. Berk and Schrieffer"^^ in- 
cluded a specific class of higher order diagrams describ- 
ing the coupling to ferromagnetic (FM) spin fluctua- 
tions, addressing almost FM metals, like Pd. They 
found that retardation is ineffective for the added dia- 
grams and that superconductivity is strongly suppressed, 
which can help to explain the cases when FM spin fluc- 
tuations are important. In another line of work going 
beyond MA the jellium model for electrons was consid- 
ered and the Eliashberg equation were solved including 
the frequency and momentum dependence of the dynam- 
ically screened Coulomb interaction. Curiously, random 
phase approximation (RPA) calculations by Rietschel 
and Sham yielded negative values for /i* when the den- 
sity parameter exceeded values > 2.5i^ The re- 
sulting unrealistic large values for Tc were interpreted as 
plasmon induced superconductivity."^^ It was shown that 
vertex corrections lead to a reduction of Tc and partly 
cure this problemi^i^^ Phenomenological models of the 
form of Kukkonen and Overhauser-^ for the screened 
Coulomb interaction mostly yield realistic positive re- 
sults for /i* ^ 0.1 and no s-wave superconductivityr^i^i 
although there also exist different conclusions.— We re- 
mark that Richardson and Ashcroft^"^ arrived at a rather 
accurate prediction for Tc in Li at ambient pressure based 
on such calculations. Nevertheless, there does not seem 
to be a conclusive study, which systematically analyzes 
higher order corrections to the unphysical RPA result. 

Rather than treating the Coulomb interaction in the 
electron gas, our approach is based on a model with local 
interactions;^ It can thus be interpreted as only consider- 
ing the second step of the MA procedure, where the pro- 
jection of the screened interaction to the phonon scale is 
considered. To be specific we take the Hubbard-Holstein 
(HH) model, where there is a local electron-electron re- 
pulsion and an electron-phonon interaction present. The 
advantage of the restriction to this model is that it can 
be analyzed by means of the dynamical mean field theory 
(DMFT)4^ Since the DMFT becomes exact in the hmit 
of large dimensions, it can provide benchmark results in- 
dependent of the interaction strength and thus allows us 
to test otherwise uncontrolled perturbative approxima- 
tions. Retardation effects encoded in the frequency de- 
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pendence of the self-energies are fully contained. DMFT 
treats both the electron-phonon and the electron-electron 
interaction non-perturbatively and therefore include con- 
tributions to /i* to all orders. However, it also includes 
renormalization effects of the phonons, the electrons and 
the electron-phonon coupling. This makes the interpre- 
tation of the results and the extraction of /i* more diffi- 
cult. In order to be able to nevertheless obtain meaning- 
ful insights, we developed a diagrammatic perturbative 
approach for the HH model in the limit of large dimen- 
sions. This allows us to see how accurate the conventional 
theory describes the benchmark results and which correc- 
tions are necessary to get a good agreement between the 
perturbative results and DMFT. We will see that higher 
order corrections to /i* enter in a modified form when 
compared to Eq. ([T]). 

The occurrence of superconductivity in the HH model 
has been analyzed theoretically beyond ME theory42rii 
Freericks and Jarrel'*" studied the suppression of the in- 
stabilities towards charge density wave (CDW) forma- 
tion and superconductivity at and away from half fill- 
ing within the framework of DMFT. Their finding of ro- 
bustness of CDW against superconductivity could be ex- 
plained within weak coupling perturbation theory with- 
out invoking corrections to the pseudopotential result in 
Eq. Functional renormalization group studies^i^ 

found spin density wave (SDW), CDW and different su- 
perconducting instabilities. The phonon scale in those 
works was however relatively large such that retardation 
effects only played a moderate role. 

In this work we show that retardation effects lead to 
the reduction /ic — > fJ-* also in the higher order calcula- 
tion, but not as efficiently as for the first order one. Non- 
perturbative DMFT calculations clarify that the per- 
turbative result is accurate up to intermediate coupling 
strength. An important conclusion is then that retarda- 
tion effects indeed lead to rather small values of /i* , even 
when contributions beyond the standard theory are con- 
sidered. For systems with sizable Coulomb interactions 
/ic, our values for fi* are however larger than in the stan- 
dard theory and therefore lead to reduced values of the 
superconducting gap and Tc- The paper is structured 
as follows: In Sec. II, the details of the model are in- 
troduced as well as some basic properties of the DMFT 
approach. In Sec. Ill, we discuss the diagrammatics for 
the calculation of Tc from the pairing equation, includ- 
ing self-energy and vertex corrections. In Sec. IV, similar 
diagrams are discussed for the calculation of the gap at 
T = 0. In Sec. V, we focus on the calculation of the 



fi* and derive analytic results including the higher order 
corrections. In Sec. VI, we put together the results from 
the perturbation theory and validate our findings up to 
intermediate couplings with the non-perturbative DMFT 
results, followed by the conclusions. 

II. MODEL AND DMFT SETUP 

The purpose of our work is to obtain generic insights 
into the behavior of a coupled electron-phonon system. 
Specifically, we employ the Hubbard-Holstein model, 

H = -^(t»j4<xCj,^ + h.c.)-l-i7^<^na (2) 

+OJ0 4b^ + gY^^b, +b^i){Y. - 1) • 

cj ^ creates an electron at lattice site i with spin cr, and 

6| a phonon with oscillator frequency wq, fii^a = c|^Ci^o-- 
The electrons interact locally with strength t/, and their 
density is coupled to an optical phonon mode with cou- 
pling constant g. The local oscillator displacement is 
related to the bosonic operators by Xi = {bi + bl)/^/2ujo, 
where h = 1, and one can define a characteristic length 
xq = l/y/ujQ for the oscillator. We have set the ionic 
mass to M — 1 in 1^. The model in Eq. Q possesses 
the minimal ingredients necessary, such as energy scales 
for electrons, phonons and their interactions. In the limit 
of large dimensions this model can be solved exactly by 
DMFT. Hence, we can provide controlled benchmark re- 
sults in this situation. The following calculations are 
based on this model in the limit of large dimension. In 
this case the self-energy is independent of the momentum 
fc, but retains the full frequency dependence. This can 
be compared with the usual application of ME theory, 
where one usually projects to the Fermi surface and only 
deals with frequency dependent quantities. 



A. Calculating Tc 

In DMFT the critical temperature can be calculated by 
analyzing the relevant susceptibility. For completeness 
we display some of the results, which we use later. The 
notation follows the one in Ref. H^. The equation for 
the s-wave superconductivity susceptibility Xq(ia;„) reads 

with Xq{i^n) = E„i,n2 Xq{i^ni , iUJui', i^n) , 



Xq{ii0niJ(^n2;i^7i) = Xq(icj„i , iw„2 ; ii:j„) -|- - Y Xq(*w„i,«i:j„3;iw„)r^PP)(ii:j„3,iw„4;icj„)xq(jw„4,ii:^„2;^ (3) 

n3,n4 

I 

r'^PP^(ia;„j , ioJn^; iwn) is the irreducible vertex in the corresponding pair propagator Xq(ia;„j ; icLi„) is 
particle-particle channel which is local in DMFT. The 

Xq(it^„i,«w„2;iw„) = 2^Gfc(«a;„JGq_fe(ia;„-iw„J5„i,„2, 
fc 

(4) 
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For special values of q and iwn we can evaluate the pair 
propagator. We are interested in the limit q — ?> and 
ibJn — > 0, and one finds 



where C(icLi„) = iLU„ + fJ- — S(iw„) and 



(5) 



de 



Pq 



CiiuJn) - e 



HT[po](C(«c^«)). (6) 



4i2 



(9) 



At half filling /i = 0, G{—iujn) = — G(iaj„), iaj„) = 
— I](ia;„). With E(iLLi„) = ia;„(l — Z(i(u;„)), where Z{iujn) 
is symmetric and real, we find 



(10) 



This expression is positive. 

We can write Eq. ([3]) as a matrix equation (omitting 
the general arguments, q,ia;„), 

x = x" + x°rx. (11) 

The instability criterion is that x diverges for some 
T, ia;„, q. This can be written as the eigenvalue equa- 
tion [x'^r]''^ = V, or with V — [x^]^^"^^, 

[xy/'r[xr^^i^i, (12) 

where the matrix is symmetric if F is symmetric. The 
relevant symmetric matrix reads, 

Mn„n, = ^ VX° ) [^^""^ i^^n, , tUJ„, ; 0)] («^«. ) , 

(13) 

and we have to find its largest eigenvalue. We have sim- 
plified the notation for the arguments for omitted the 
q = 0-label. In a DMFT calculation for a given temper- 
ature T, we first determine S(«aj„) and the full particle- 
particle irreducible vertex F^PP-'(jLij„j, iw^j; 0). The pair 
propagator x'^(iaj„j) is obtained from Eq. ([S]). Then we 
can search for the largest eigenvalue of the symmetric 
matrix in Eq. (jl3p for the instability criterion. 



B. Calculations in the superconducting phase 

We also perform calculations in the superconducting 
phase. We work in Nambu space with matrices then. 
The local lattice Green's functions have the form 



Gn(iL.„) = AGHT[po](e+) + BGin[po]{e-) (14) 



and 



One has for the semi-elliptic DOS 

D 

mPoKz) = Jde^^ ^(z-sgn(Im(z))v/z2-4t2), 

-D 

(7) 

where the square root of a complex number w is given by 
y^e*'^/^, where ip = [0, 27r), such that the imaginary part 
of ^/w is positive. At half filling Giiujn) and S(iaj„) are 
purely imaginary functions. The phonon Green's func- 
tion D{iuj.ra), 

D{iujyny^ = D°{iuJrny^ - Sph(ia;„i), (8) 

where D'^(iujm) — 2ajo/[(it^m)^ ~ ^q], and its self-energy 
I]ph(ia;m) are real functions. In the non-interacting case 
the Green's function reads. 



G2i{iL0n) = Aj.HT[po](e+) + Bf^HT[po](£- 



(15) 



with Ag = (C2(iw„) + e+{iujn))/{e+{iujn) - £-(iw„)), 
Bg = {C,2{ii^n) + e-{iujn))/{e-{iLOn) - e+{iujn)) 
Af = i;2i(«a;„)/(£+(iw„) - £_(ia;„)), and Bp = 
S2i(iw„)/(£-(ia;„) - e+(iw„)), where 

Ci(ia;„) - (2{ii^n) , 
= 2 ^ 

2'\/(Ci(«'^") + C2(«a;„))2 _ 4S2i(iw„)Si2(iw„), 

withCi(z) = Eii(z) andC2(z) = z — ^j,—T,22{z)- We 

have Gi2(«w„) = G2i(iw„) and G22(«Wn) = -Gii(-ia;„) 
for the Nambu Green's functions. We use Ei2(ia;„) — 
T,2i{iujn) and S]22(^w„) = — I]ii(— ja;„) for the self- 
energies. This can be deduced from the properties of the 
corresponding Green's functions including the assump- 
tion of time-reversal symmetry. At half filling Gii(iaj„) 
and I]ii(ia;„) are imaginary functions, whereas G2i(ii^„) 
and Yi2i{iuJn) and D{iujm) and "Epi^^iujm) are real func- 
tions. 

In the NRG approach we calculate the self-energy 
matrix for the effective impurity model from the ma- 
trix of higher Green's functions F_{uj) with Fii{uj) — 

((cd,t'^i;4,t))"' ■^12('^) = {{cd,tni;cd.i))u,, F2i{uj) = 
-((ci;'^t;4,t))" and F22(cj) = Crf,;))^. For 

the phonon part we use Mii(c<;) = {{cd.-\-{b + b^); c^^ ^))^, 
Mi2{uj) = {{cdA{h + h^)-Cd,i))u., M2i{uj) = -{{c\^{h + 

b'^)\c\^))u, and M22(^^) = -ic\ Ab + h^)-Cd.x))u- In the 
NRG we calculate Mn and M21 directly on the real 
axis from the Lehman spectral representation. The oth- 
ers follow from Mi2(a;) = — M2i(— w)* and M22(w) = 
Mii(— cj)*. We can define the self-energy matrix by 



E(w) = UF{uj)G{uj)-^ + gM{u)G{uj)- 



(16) 



For self-consistency the local lattice Green's functions 
Giuj) has to be equal to the impurity Green's function, 
Gioj) — G{u}), where 



g-\uJ)=LJl2 + ^lT3-K{io)^U^), 



(17) 



with the matrix K_[ijj) describing the effective medium. 
We can take the form of the effective impurity model 
to correspond to an Anderson- Holstein impurity model— 
and calculations are carried out as detailed, for instance 
in Ref. 47. We solve the effective impurity problem with 



the numerical renormalization group"*^'^^ (NRG) adapted 
to the case with symmetry breaking. The NRG has been 
shown to be very successful for calculating the local dy- 
namic response functions, and we use the approach^-2i^ 
based on complete basis set proposed by Anders and 
Schillerj^ For the logarithmic discretization parameter 
we take the value A = 1.8 and keep about 1000 states 
at each iteration. The initial bosonic Hilbert space is re- 
stricted to a maximum of 50 states. We will mainly con- 
sider two cases: (i) constant density of states po = 1/14^, 
where W is the bandwidth and (ii) the semi-elliptic DOS 
poie) = VW~e^/{2TTt^) with W ^ 4t ^ 2D. 

III. DIAGRAMMATIC CALCULATION FOR n 

In this section we first show how the standard approach 
to conventional superconductivity, ME theory, would be 
applied to the model under consideration, and which di- 
agrams are included. To determine we use the in- 
stability criterion, Eq. (I12p . which is equivalent to the 
linearized ME equations. 

A. Standard diagrammatics, ME theory 

In ME theory the irreducible vertex in Eq. is given 
by the full phonon propagator 

r'^PP)(ia;„i,«i:j„2;iw„ = 0) = -g^D{iujn^ ~ ii^n^)- (18) 

Diagrammatically this is depicted in Fig. [1] (a). 




(a) (b) 



FIG. 1: Contributions to the irreducible vertex (a) electron- 
phonon part, (b) first order in U . 

Due to Migdal's theorem^ ^ other vertex corrections in- 
cluding the phonon propagator are neglected. In gen- 
eral, the ME theory is a self-consistent calculation, and 
the electronic self-energy reads [see Fig. [5] (a)], 

S(iw„) = ^ G(iw,„ -I- iuJri)D{iuj„i)- (19) 

m 

The phonon self-energy reads [see Fig. [5] (b)], 

5^ph(iw,„) = G{iujn)G{iuj,n + iuJn). (20) 
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D 




(a) (b) 



FIG. 2: Diagrams for (a) electronic self-energy and (b) the 
phonon self-energy. 

We define as the relevant renormalized phonon scale, 
which is extracted from the peak position of the full 
phonon spectral function p^{u!). The electronic Green's 
function G{iuin) is determined as in Eq. ([5]). An appro- 
priate definition of the coupling constant A is 

oo 

A = 2|dc.^^, (21) 



where we the pairing function is defined via a Fermi sur- 
face average, 

a^F{uj) ^\Y1 <k'F{uj)S{ek - p)d{ek' - m) (22) 
with 

al^,F{oj)^po\gk,k'\^ptk'i^)- (23) 

For the Holstein model using the free spectral function 
we have a^i^(ti;) = pog'^[6{uj — uo) — 6{uj + luq)], such that 
we obtain Aq — po2g'^/uio, purely in terms of bare param- 
eters. However, more relevant is the effective coupling A 
for the Holstein model is given by the full renormalized 
phonon propagator, 

X^2pog^ fdu ^-^^-pog^D{Q). (24) 



The dimensionless quantity for the Coulomb interaction 
corresponding to Aq is pc = PqU. 

In a self-consistent numerical calculation for a given 
temperature T, we first determine S(iaj„) and Sph(itLi,„) 
by iterating Eqs. and ([^0]) and thus G{iuJn) and 

D{iuJm). Then we determine from Eq. 

D{iujm) is used in Eq. to determine the irreducible 
vertex. Then we calculate the largest eigenvalue of the 
symmetric matrix in Eq. (jl3p for the instability criterion. 
Instead of calculating the phonon Green's function self- 
consistently we can also take it as an input. Using the 
DMFT results for D{iujm) we found in previous worbi^ 
that Tc obtained from this procedure agrees well with 
the full DMFT result as long as the renormalized Migdal 
condition, Xujq/W small, is satisfied. In the DMFT cal- 
culation r'PP' and S(iw„) contain all higher order cor- 
rections. 

In the usual theory the Coulomb repulsion is included 
directly up to first order or as an additional empirical 
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parameter ^* . The vertex has then two contributions, 
one as before due to the electron-phonon interaction, and 
the second one from the Coulomb repulsion [see Fig. [1] 
(b)l, 

r<^PP^(ia;„i, iwn^; iwn = 0) = -g'^D[iuJm - ^^nj - U. 

(25) 

The eigenvalue equation (|12p with the pairing vertex 
(|25p can be solved analytically with approximations sim- 
ilar to the ones by McMillan.^^ This yields 



2Cw?; 



■ exp 



A 



(26) 



with C — e^ PS 1.78 with the Euler-Mascheroni con- 
stant 7 = 0.57721. We have introduced Z = 1 — 
limtj^o 'S{iio)/iuj and /i* is given by Eq. ([T]) with E^i = D 



and 



■^ph 



McMillan^ or Allen and Dynes^'', 



This corresponds to the result by 



(Uj) 1.0i(l+\) 
2 i_p A-M* (1 + 0.62A) 

1.2 



(27) 



where Z = 1 + A is used. The essential feature is that 
the Coulomb repulsion, which is the same for all ia;„ 
is effectively reduced from /ij. to /z* when counteract- 
ing the electron-phonon attraction with strength A. This 
shows how retardation effects assist the electron-phonon 
induced superconductivity by suppressing the detrimen- 
tal effects due to the Coulomb repulsion. 

In summary, the diagrams in Figs. [T]and[2]are the ones 
included in the standard theory of superconductivity. 
Usually, the phonons are not calculated self-consistently 
but rather taken as an input, for instance from DFT cal- 
culations or from experiment. Also /i* is usually not cal- 
culated but used as a fitting parameters. In the following 
we consider higher order corrections to the standard ap- 
proach. 



(b) Contributions r|j^'' which include the bare 
vertex Fg^P^ and higher order terms purely due 
to U. 

(c) Higher order corrections T^^^ due to mixed 



(cp) 



terms of g and U, not including Fq 

2. Contributions to the phonon self-energy Sph: 

(a) Higher order contributions to Sph, which can 

be written in terms of Fg''^-', i.e., purely due 
to g. 

(b) Contributions to Eph due to mixed terms of 
g and U. These can be expressed in terms of 

3. Contributions to the electron self-energy E: 

(a) Contributions to S, which can be written in 
terms of Fg°'^', i.e. purely due to g. 

(b) Contributions to E purely due to U. 

(c) Contributions to E due to mixed terms, which 
can be expressed in terms of F^^^'' or T^g^- 

(d) Contributions to E due to mixed terms, which 
can not be written in terms of F^^p^. 

4. Contributions to the irreducible vertex F^pp^: 

(a) Contributions to F*^pp) purely due to U. 

(b) Contributions to F^pp^ due to g and U, which 
can be written in terms of full propagators and 
the electron-phonon vertex F'^°p\ 

(c) Contributions to F^pp-' due to g and U, which 
can not be written in terms of F*^°p). 



B. Higher order terms 

The perturbation expansion for two different interac- 
tions is rather involved, since terms of each perturbation 
series as well as mixtures can appear. In a skeleton ex- 
pansion higher order contributions can be grouped into 
the following terms for self-energies and vertex functions: 

1. Contributions to the full electron-phonon vertex 



p(op) _ r(°p) 

IS F^ g: 



F 



(op) 



n(op) 



where the bare vertex 



(a) Higher order corrections Fg°^'' purely due to g 
not including Fq'^''^. 



We assume in the following that the parameters are cho- 
sen such that higher order contributions of the type Fg''^'' , 
i.e., 1(a), 2(a) and 3(a), are small due to Migdal's the- 
orem. By taking the full phonon propagator from the 
DMFT calculation as in previous workfii we avoid con- 
sidering in detail effects of 2(b), and rather assume that 
we can include the correct phonon propagator. We focus 
on contributions 1(b), 3(b,c), and 4(a)(b), in the follow- 
ing, which are the main contributions in the low order 
perturbation theory. Contributions to the type 1(c), 3(d) 
and 4(c) include higher order diagrams and will not be 
considered explicitely here. 

We can generally write contributions to the pairing 
vertex of the form 4(b) as 



r(pp)(i 



,IUJ„ 



0) 



^(cp; 



,lUJr 



,)D{iujn 



-IUJ„ 



(28) 
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where T^j^^^ includes g and all corrections due to U, see 
Fig.H 



FIG. 3: Contributions to the electron-phonon vertex r[°^\ 

For r[;P^(zt^ iwnj, iwni is the ingoing electronic fre- 
quency, iw„2 the outgoing one and the bosonic one is 
iujn-^ — iuJn^. D(iujn) is the full propagator including cor- 
rections due to U and g. Diagrammatically, these contri- 
butions to r^PP'' are displayed in Fig. 5] (a). 



1^,11 1 





(a) 



FIG. 4: (a) Higher order terms involving vertex corrections 
of the electron-phonon vertex, (b) Electronic self-energy in- 
cluding vertex corrections for the electron-phonon vertex. 



The higher order corrections included in Eq. (|28p can 
be seen as a redefinition of A defined in (|24)) . The ver- 
tex r^^^-* does not vary much up to the small phonon 
scale, such that we can write approximately with = 



A~2po[5l' /dw 



(29) 



where p^{uj) can also include self-energy corrections due 
to U. Such a redefined A enters the approximate results 
for Tc such as Eq. (pS)) . We can also take into account the 

effect that F^^^'' is a function of frequency and average 
over a range given by the phonon Green's function such 
that 

1 



[9 



r-\2 



Nd 



f[;p)(0,zl.„,)2^(zl.„J, (30) 



where Nd = J2n2 ^(«^«2)- 

We now consider the diagrams contributing to the 

electron-phonon vertex T^^K Some of them are shown 
in Fig. [3l To order Ug we have 

r[jP'^''(ia;„^,iw„2) = g[l + UIl{iuj„^ - zWnJ], (31) 
where the particle-hole bubble is given by 
1 



iujm)G{iuj„i)- 



(32) 



n(0) can be evaluated analytically at T = and one finds 
for the constant density of states 11(0) = — 21og(2)po and 
for the semi-elliptic DOS 11(0) = ^|po- It turns out that 
a good empirical form for n(ia;) is given by 



-apo 



l + bi\oj\ + b2UJ^' 



(33) 



If G{iujn) is the non-interacting Green's function in Eq. 
([5^ then a = 4/3 for the semi-elliptic DOS. In this case 
we find that bi — 1.3/D and 62 = 1.63/13^ give a good 
fit. 

We can sum up a whole series of diagrams of the type 
discussed in Eq. (PT|) which corresponds to screening on 
the RPA level. 



T^(cp,RPA)/. 
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1 - c/n(zw„, 



(34) 



Note the absence of a factor 2 in the denominator for the 
Hubbard interaction. Since 11 < these diagrams lead to 
an effective reduction of the electron-phonon coupling. If 
one assumes that 11 does not change much on the phonon 
scale enforced via D{iujjyi), then one could approximate 
g'' = g/[l-;7n(0)], and for the Bethe lattice 5'' = g/[l + 
4/^c/3]. 

There are three diagrams to order in addition to 
the screening term to correct the vertex (see the third 
diagram in the Fig. [3] and Fig. [5|). 




(b) 



(c) 



FIG. 5: Vertex corrections to the electron-phonon vertex to 
order gU^. 

The first term. Fig. [5] (a), reads, 

^ k 

xG{iujk + iujm). (35) 
The second one. Fig. [5](b), is of the form 

p(op,2,2)^^^^^^^^^ = Upp{iuJk)G{iuJk - iuJn) X 

^ k 

xG{iujk - iujm), (36) 
where we have introduced the particle-particle bubble, 

Upp{iUn) = ^ X! ^(^'^'i ^ iu!„i)G{iuj,n). (37) 
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At half filling we have G{—iLJ,n) = —G{iujjn), which im- 
plies Iipp{iujn) = ~n(iw„). With this one obtains 



We plot results for {cf / g)^ in Fig. [SI 



U 



{iuJn,iUJm.) = -r 



(op, 24) 
U 



(38) 



and the first diagram is canceled. The third diagram, 
Fig. [5] (c), is like the first one r[j'''^''^''(iaj„, iwm) = 

r[j^'^'^''(jw„, icjm)- So altogether we have just one con- 
tribution of the form ([55]) at half filling. Diagrams of 
these type together with the RPA screening series were 
discussed by Huang et al.^^ for the two dimensional 
Hubbard-Holstein model. The full electron-phonon ver- 
tex was calculated with QMC. It was found that for val- 
ues of U up to half the bandwidth the diagrammatics 
and the non-perturbative result are in good agreement. 
However, for large values of U , the perturbation theory 
breaks down. 

As can be seen from this analysis of the vertex correc- 
tion, there is a considerable effect to suppress the cou- 
pling g when U is finite. Including diagrams in Eq. (j34p 
and p5p we can write the approximate form. 



g' 



1 



where a = 



g 1 + afic 
-H(0)/po and 



(39) 



(40) 



In the simplest case for free Green's functions and a semi- 
elliptic DOS we have a = 4/3 and b — 0.8237 numerically. 
In a calculation strictly up to second order in U we have 
instead of Eq. 



'- = 1 




afic + (a^ - b)^j.l 



(41) 



FIG. 6: (Color online) {g^ / g)"^ as a function of /ic, full line 
according to Eq. (|39|l and dashed line according to Eq. (|4ip . 

We see that there occurs a substantial suppression, such 
that already for /Zc = 0.25 the quantity [g^ / gY is reduced 
to about a half of its value for jic — 0. Since a? > b, 
the second term in Eq. (|411) which dominates for larger 
/ic, leads to an upturn of the result. The result in Eq. 
([5^ overestimate the suppression effect for large values of 
fic!^ For the interacting system, the coefficients a and b 
tend to be smaller than the values in the non-interacting 
limit. 



The vertex corrections T^^"^ also enter the electronic 
self-energy as one class of diagrams contributing to 3(c). 
This can be written as 



S(ia;„) 



iuJrn + i^n)G{iuJm + iuJn)D{iuJri 



+ IbJn, l(^n 



(42) 



This is shown diagrammatically in Fig. [4] (b). As dis- 
cussed above, there are also other types of mixed dia- 
grams 3(d), which cannot be written in the form of Eq. 

621). 



C. Higher order corrections from purely U 

We now deal with the higher order corrections purely 
from U, i.e., of type 4(a). We will restrict our attention 
only to the terms second order in U. The correspond- 
ing contribution (crossed diagram) to the pairing vertex 



reads 

r(PP)(iw„,,iw„,;ia;„ = 0) = C/2H(za;„, za;„J. (43) 

The diagram is depicted in Fig. [7] (a). 
A naive way of taking this into consideration would be to 
assume that H(iaj„^ -t-iwn^) varies little for the scales un- 
der considerations, such that for all frequencies we can as- 
sume — poC^^II(O) = ais^. Then this term can be treated 
in the same way as the first order term in U, where we 
can simply write ^ fie = IJ-c + '^IJ'c- This quantity is 
then subject to the full retardation effects and becomes 
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(a) (b) 



D 




(a) (b) 



FIG. 8; Diagrams for the off-diagonal self-energy (a) from 
electron-phonon interaction and (b) for the Coulomb repul- 
sion with the notation F = G21. 



FIG. 7: (a) Higher order diagram vertex contribution from 
the Coulomb repulsion, (b) Second order in U diagram for 
electronic self-energy. 



(44) 



However, as analyzed in detail in Sec. |V] retardation ef- 
fects are less efficient for the second order term when the 
decay of n(iw) is taken into account properly. 

For the self-energy contributions of type 3(b), we will 
only consider the standard second order diagram de- 
picted in Fig. [71(b), 



El 



"1 ,"2 



G(ia;„-ia;„JG'(itJ„j+iw„JG'(iw„J. 

(45) 



IV. DIAGRAMMATIC CALCULATION FOR 
THE SUPERCONDUCTING STATE 

In this section we consider calculations in the super- 
conducting state. We first present the application of 
the standard theory to the HH model and then discuss 
higher order corrections. These calculations allow us, 
for instance, to study the gap at T = 0. The relevant 
quantities are matrices in Nambu space. We have to cal- 
culate the diagonal and off-diagonal components of the 
self-energy. First we present the diagrams usually used 
in the standard ME approach. 



This is depicted in Fig. [5] (a) with the off-diagonal self- 
energy marked as F. 

In the limit of low temperature, T — 0, we use 



00 

1 . 1 f 



(48) 



For the off-diagonal self-energy [see Fig. [8] (b)] we have 
the Coulomb contribution (77 — !■ 0) 



S2i(ia;„) 



^^G2i(zc.™)e'""''. 



(49) 



Since G2i{i'-^m) decays sufficiently rapidly the factor 
e'"™'' is usually not needed. 

These are the equations, which are taken into account 
in the standard theory. Often, the Coulomb term is then 
projected onto the phonon scale (Coulomb pseudopoten- 
tial), where only the reduced value /i* enters. Making a 
few simplifying assumption (see Sec. IV B[) we can derive 
for the spectral gap at T = 0, Agp ~ Y.2i{0)/Z, 



2ujQe 



(50) 



This result is very similar to the one for Tc- Equations 
(|47p. and (P^ correspond to the standard approach 
for the theory of conventional superconductivity^^!^ 



A. Standard diagrammatics, ME theory 

To lowest order the diagonal self-energy is given by 
Eii(icj„) = -— ^ Gii{iujjn + it^7i)D{iu^), (46) 

m 

and the off-diagonal self-energy reads. 



S2l(«W„) = ^ X! '^2l(iWm + ii^n)D{iuJm)- (47) 



B. Higher order corrections 

In a similar way as for the analysis of the pairing vertex 
in Sec. IIIH we can classify terms for the self-energy into 
different contributions. We will follow the same logic as 
in the calculation for and focus on the same type of 
diagrams as before. Here we consider a subclass of type 
3(c), which includes Coulomb vertex corrections of the 
electron-phonon vertex of the type T^^^\ It reads. 
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Sii(ia;„) = -^'^T'j^^\iuJn,iujjn + iw„)Gii(iw™ + ja;„)L'(ia;™)r[jP^ (iw™ + ia;„,za;„). (51) 



Assuming a slow variation of Tu on the phonon, scale a good approximation is of the form, 



(52) 



Similarly we have for the off-diagonal self-energy, 

S2i(«w„) = i^r|^P^(ia;„,iw,„ -f iw„)G2i (iwm -I- iw„)£)(iw„i)r[jP^(ia;,„ + zw„,iw„). 



(53) 



I 

For {itOn, i^m) we consider the same diagrams as in propagator and hence it contains corrections due to U as 
Sec. mil B. We assume that D{iuJm) is the full phonon well. 

I 



C. Higher order corrections from purely U 



Of the contributions of the type 3(b) we discuss all terms to second order in U. In the Nambu perturbation theory 
ones has (cf. Ref. i56. ). 



(54) 



and 



Ylf^^^iiujn) ^ ^ G'2i(ia;„ + iw„JGi2(jwm2)G22(jwmi + ia;,„J. 



(55) 



For the off-diagonal part we have 

^2l',C/(*'^") = ~(^) X! G2l(ia;„ + itJ™JGi2(iWm2)Gi2(iWmi 



(56) 



v^(2,2)/. N 



^ Gii(za;„ -I- iwmJG22(jwmJGi2(ia;mi + iwmj. 



(57) 



r 



The first diagram E|^^']^(icj„) is the well known U^- 
diagram which gives the first dynamic correction, also 
in the normal phase as in Eq. (|45p . In comparison 

f2 2) 

jj{iujn) gives a smaller contribution as it is propor- 

f 2 1) 

tional to two off-diagonal Green's function. '[/(^'^n) 

(2 2) 

is also comparably small, but Yl^i ui'^'^n) gives a sizeable 
reduction of the superconducting state. We can see this 
by writing it as 



^2l',t/(*'^n) = -t/^^y^n(za;n+i^mi)G2l(i^mi). (58) 



where n(i(u;„) is given in Eq. (|32p with G{iujm) 

Gii{iu}m)- 

For small iw„ a crude approximation is to write 



^2i$i0) ^ -C/'cin(0)i^G2i(*c.™J, (59) 



where it is assumed that G2i(iw„ij) is only finite in small 
interval such that we can take n(jcj„) constant. Hence 
we find a direct correction to the term 'S^^jj. With the 
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result for 11(0), one has then approximately 



4i,u + S^"c^(0) ^ C/(l + oMc)^ X. G2i(*wmj, (60) 

mi 

This is the same effect that was discussed for and 
the crossed diagram [see Fig. [7] (a)] where naively the 
effective fic becomes flc — /^c(l+iMc)- In a more accurate 
treatment, we have to take the frequency dependence into 
account, and we will see that this leads to modifications. 



.(2,2) 
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V. ANALYTIC AND NUMERICAL RESULTS 

FOR 11* 

In the last two sections we have analyzed the diagram- 
matic expansion for the model in Eq. ^ and discussed 
certain types of diagrams. In this section we want to 
specifically study the pseudopotential effect without in- 
cluding all other corrections. We are interested whether 
the first and second order calculations give qualitatively 
different results for /x*. We present a combination of 
analytical and numerical arguments. In the literature, 
there exist a number of ways to calculate fi* for a given 
microscopic model. One early approach by Bogoliubov 
et al. is based on an integral equation for the Coulomb 
part in the pairing equation, Morel and Andersoi*ii 
gave an approximate solution for the Migdal-Eliashberg 
equations including a screened Coulomb repulsion in the 
T = formalism. The pseudopotential /i* also appears 
naturally when superconducting pairing instabilities are 
studied in a renormalization group framework4^iS In the 
following we first calculate fj,* directly by projecting the 
pairing matrix to low energy. Analytic and numerical 
results are compared. Then in a second approach we 
calculate Agp from an approximate solution of the self- 
consistency equation and thus extract fi* . Finally we also 
compute numerical results for Tc obtained from the in- 
stability condition and analyze these results in terms of 



A. Projection scheme 

The starting point for the projection approach is the 
pairing matrix in the form 

^nm — ^nm ^^nm^ (^I) 



where Mnm is given in Eq. (|T3|) . This matrix becomes 
singular at Tc. We introduce the "low-energy part" 



Alow _ A 



'^nl-l'^ml>'^pli 



n, m such that |i^n|, I'^ml < '^ph, and A is the matrix that 
is left after the blocks given by |a;„| < Wph or |aj,„| < 
Wph were removed. If is not singular, Eq. (|6ip and 
Eq. (j62p become singular for the same parameters. This 



way we can reduce the matrix size so that it only includes 
frequencies for which the electron-phonon interaction is 
important. The "folding in" of larger frequencies then 
describes how retardation effects reduce the effects of the 
Coulomb repulsion on low frequency properties. 

We first consider the lowest order term of F'pp) in U 
namely F^pp^ = —U. We want to focus on the dependence 
on the half-band width D and assume that the density 
of states scales as 



Pie) 



(63) 



where p is independent of D. For simplicity we assume 
in the following that p{x) = 1 is a constant for \x\ < 1/2 
and zero otherwise, such that po = 1/(21?). It is then a 
rather good approximation to write 



poTr/lujnl, if |a;„| < D; 
0, otherwise. 



The matrix A then takes the form 

TT 



(64) 



(65) 



/3y/\uJnUJra\ 

which is separable and can be inverted exactly. We obtain 



[A ]nm ^nm 



(66) 

where involves a summation over |ajfc| > Wph. Re- 
placing the summations in Eqs. (|62[ I66p by integrals, we 
find 



4I0W _ c I 



l3\/\uJnUJ„ 



;3^|w„a;„| 1 + Aiclog(£'/wph) 



(67) 



Comparison of Eq. (|65|) and (j67| leads to the Coulomb 
pseudopotential, fj.c ^ p*, as given in in Eq. Ilj, and it 
is the result of Morel and Anderson. 

We next consider the second order term of F^pp^ in U 
[Eq. (|43p]. Because of the form of Il{iujn + iuj„i) it is then 
not possible to invert A in Eq. ([S^ analytically. For the 
density of states [Eq. we can rewrite n(ia;„ -|- iujm) 
as 



(68) 



where / is independent of D. As in Eq. ([33|) it is quite 
accurate to approximate / as 



/(•^) = 



1 + 61x1 + cx^ ' 



(69) 



where a = 1.38, 6 — 2 and c — 5 are suitable values for 
the constant DOS. 
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We now want to calculate A^™_ in Eq. to third 
order in U. We use the result for in Eq. which 
neglects the second order term and is therefore only cor- 
rect to order U. However, since the off-diagonal terms of 
A in Eq. are of the order U, the final result is of the 
order U^. We obtain 



nm 



Aln[A 



nniAmi 



2D J 



'i-Y- 

\I3J 1 



(70) 



where we focus on the diagonal result for the lowest fre- 
quency LUi. Using the definition of x'^{iu!n) in Eq. (O, 
not the approximation in Eq. (I64[) , and applying the limit 
i^ph <^ D, we calculate the sums 



E 



<i'/5 /DAi, 
log 

2D ^Vujph 



(71) 



where Aq « 1.00, Ai « 0.32 and A2 « 0.20. For /c = 1 
and 2 the denominator of / reduces the integral for large 
\u)n\, which effectively reduces the band width by a factor 
Afe. This reduction is naturally larger for /c = 2 than for 
fc = 1. 

We now make an ansatz for /i* 



22 



(72) 



along the lines of the Morel-Anderson form, but includ- 
ing a second order term and a corresponding logarithm. 
Based on Eq. ([7T|) we expect the effective band width 
to be smaller for the second order term and we there- 
fore allow for a different multiplying factor A22 in the 
logarithm. Since the result in Eq. ([70)) is correct to 
order the factor A22 in Eq. ([7^ can be identified, 
A22 = Aj/Ao « 0.10. The ansatz in Eq. is then also 
correct to order U^. 

Fig. [0] shows the results obtained by performing the 
calculations in Eq. (p^ numerically using the first or 
second order result in U for F^pp^ and without intro- 
ducing the approximation in Eq. (j64p . The figure also 
shows the analytical result in Eq. ([7^ . The second or- 
der result is clearly larger than the first order result. For 
/ic ^ 0.5, Eq. (I72p describes the full second order calcula- 
tion rather well, while for larger fic there are corrections 
to the analytic result which make fi* still larger compared 
to Eq. ([72)1 . 

It is interesting to discuss the origin of the factor A22 
in Eq. ([7^ . In the present language the term log(L'/ajph) 
in the first order calculation arises from "folding in" large 



frequency contributions from U in Eq. (j62p . which ex- 
tend to approximately u ^ D. The second order term 
in U is "folded in" in a similar way. However, as shown 
by Eq. (j69p , the contribution from large frequencies is re- 
duced, leading to the effective band width being reduced 
by some factor for the second order contribution. This 
reduction factor is surprisingly large (1/10). The origin 
can be seen in Eq. (|70l) . where the relevant terms are 
a product of the first and second order contributions in 
the first term, containing a factor two, and the first order 
contribution in the second term. The prefactor two of the 
logarithm in the first contribution leads to the factor A\, 
and thereby a very small factor. 

This shows that the second order terms contribute to 
retardation effects, like the first order contribution. How- 
ever, the reduction is substantially less efficient, as de- 
scribed by the small factor ^22- The band width there- 
fore has to be very large to make this contribution to 
retardation efficient. 



Second Order. Calc. 
Second Order. Analyt, 
First order. Calc. 




FIG. 9: The pseudopotential /i* as a function of fic for 
D/ujph ~ 100 and fiujph ~ 240. The figure shows the cal- 
culated results using both the first order and first plus second 
order result for V^'^'^' as well as the approximation in Eq. (|72|l . 

It is interesting to consider the second order contribution 
alone. We can then calculate results accurate up to order 
C/^. The matrix A^^ can be approximated by a unit 
matrix. Making an ansatz for fi* and identifying with 
the fourth order result, we obtain 



(73) 



where A2 = 0.20 was given above. The result of a full cal- 
culation is compared with the analytical approximation 
[Eq. (|73p [ in Fig.[Tni Also in this case the analytical result 
agrees rather well with the full calculation for /ic ^ 0.5. 
Finally, we compare the first and second order result as 
a function of log(£'/wph) in Fig. [TT] Different values for 
fic were used in the first and second order calculations, 
and chosen such that both calculations gave the same 
fj,* for D/ojph = 10. It is quite interesting that ^* then 
decreases in a very similar way in the two cases. 
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FIG. 10: (Color online) function of fic for D/ujph = 100 

and l3ujph ~ 240. The figure shows the calculated value using 
just the second order result for r'^^' as well as the analytical 
approximation in Eq. (|73p . 




0.13 I • ■ • • ■ • 1 

2 2.5 3 3.5 4 4.5 5 5.5 

log(D/o)pf,) 

FIG. 11: (Color online) function of log(i3/a;ph) for 

/3ajph — 240 according to a first {fici ~ 0.5) and second (^c2 = 
0.2756) order calculation, /ici and fic2 were chosen so that the 
same ^"^ was obtained in the two calculations for D/uph = 10. 



We are now in the position to discuss the results for 
We have performed calculations using the first or sec- 
ond order expression for F^pp'. We then either calcu- 
lated Eq. (|62p numerically, without any further approx- 
imations. Alternatively, we performed analytical calcu- 
lations involving the neglect of additional higher order 
terms in U. 

When the first order result for F^pp' is used, the ana- 
lytical calculations can be performed without further ap- 
proximations. This leads to the Morel-Anderson result, 
which has two important features. When D/ojph oo, 
^* — > 0, and when [/ — oo, ^* — > l/log(£>/a;ph) stays 
finite and saturates. Is this still true when the second 
order contribution to F^pp^ is included? 

The analytical results in Eqs. (|72l73p have these prop- 
erties. However, these results were not derived but are 
ansatze inspired by the Morel-Anderson result and ad- 
justed so that they agree with the analytical calculations 
to low order in U. Actually, Figs. 191 liH show that although 



retardation effects strongly reduce n*, there is no sign 
that the values saturate as U becomes very large. In this 
sense there is an important difference from the Morel- 
Anderson result. For these large values of U, higher or- 
der effects in F'^pp^ become important, and it is not clear 
how these influence the conclusions. 

The second issue is how fi* is influenced when D/ojph 
becomes very large. As is clear from Eqs. (|72l73p . retar- 
dation effects also reduce the second order contribution 
to /i*. This is also seen in Fig.[TT] However, the effective 
band width is smaller due to the frequency dependence of 
the second order contribution and retardation effects are 
less efficient. Nevertheless, the second order contribution 
goes to zero as D/ujp\^ — >■ oo. In this context Fig. Illl mav 
seem surprising. One might have expected the second 
order contribution to drop more slowly with D/ujp\^. To 
understand this one can study Eqs. (I72I73P and find the 
/ici and /ic2 which lead to the same /i* in the first and 
second order calculation. Because of the less efficient re- 
tardation effects for the second order term, /ic2 has to be 
chosen smaller than would otherwise have been the case. 
The criteria for the choice of /Xc2 is, however, independent 
of D. Thus the two curves in Fig. [TT] should be identi- 
cal according to Eqs. (|72l73p . The small deviation is due 
to the inaccuracies of these analytical results for finite 
U . Fig. [TI] nevertheless, nicely illustrates how both the 
first and second order contributions are systematically 
reduced as D/cUph is increased. 

B. Calculations for the gap at T = 

In this section we conduct a complimentary analysis 
to extract results for /i*. We calculate the spectral gap 
of the superconductor Agp at T = 0, and show how /.t* 
enters naturally in the analytical description. We present 
this analysis to order [/^. The approach here is similar 
to the original work by Morel and Anderson, which in- 
cluded only the first order term in U and was carried out 
in the T = formalism. Here, we work on the imag- 
inary axis in the limit T — > 0. Starting point is the 
self-consistency equation for the off-diagonal self-energy, 

S2i(ii^„) = ^'^G2i{iuj„i)K{iujn,ii^m), (74) 

m 

where the kernel K{iujn,iuJm) includes the following 
terms, 

K{iujn, iujm) ~ g^D{iujn — ii^m) + U — U^Il{iuJn + i^m)- 

(75) 

The second and the third term are as given in Eq. p9)) 
and Eq. (|58p . The expression for G2i{i(^m) was given in 
Eq. and a semi-elliptic DOS is used in the following. 
We do not consider vertex corrections of the electron- 
phonon vertex here, and we simply use the form, 

g'Di^cu„) ^ i— ^. (76) 
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The effect of the diagonal self-energy is taken into ac- 
count in the analytical calculations for completeness with 
a Z-factor for small frequencies, |a;„| < Wph. In the nu- 
merical calculations in this section it is neglected. 

The self-consistency equation ([71)) can be solved nu- 
merically by iteration to find a solution for S]2i(iw„). For 
an analytical solution, we need to make some approxima- 
tions. At half filling we use for the Green's function for 

l^^nl < Wph, 



6*21 («W,0 ~ - 

for Wph < |a;„| < D, 



1 



G2l(jw„) 



1 S2i(w„) 
t 



and for 



> D 



G2l{iUJn) 



S2i(it^„) 



0. 



(77) 



(78) 



(79) 



As discussed, Il(iuj) is well approximated by the form 
given in Eq. (|33|) . As noted by Morel and Andersortii the 
w-dependence of the off-diagonal self-energy is very sim- 
ilar to the one of the pairing kernel K{iujn, itOm)- Hence, 
a suitable ansatz for the off-diagonal self-energy is, 

A2 Ai - A2 - A3 



S2l(iw) = A3 + 



1 + bi\uj\ + b2Uj^ 



1 



(80) 



It contains the parameters A^. We assume that 61 and 
&2 take the same values as what is found for Il{iuj) in 
Eq. By comparing the w-dependence with the full 

numerical solution we find reasonable accuracy for this 
assumption. In general we have to solve then for the 
parameters three Ai, A2, and A3 by evaluating the self- 
consistency equation at suitable values of iu. Unfortu- 
nately, the general case is algebraically very involved. It 
will be discussed to some detail in the appendix. Here 
we only treat the first and purely second order cases to 
see the major effects. 

For the first order case, we set A2 = and omit the 
C/^-term in Eq. ([751) . We use the self-consistency equa- 
tions E2i(0) = Ai, and T,2i{iD) ~ A3. When evaluat- 
ing the S2i(0) and T,2i{iD) according to Eq. ([7i)) . we 
approximate the integrals in order to find an analytic so- 
lution (see appendix). We also assume A,; ^ cjph ^ D 
for simplification. Then for a self-consistent solution the 
parameters Ai and A3 have to satisfy the equations 



Ai 



A 



Ailog 



ph 



A3 - -n 



Z -^-""V Ai 
Ai /2Za;ph 



aoAA3 - ^cA3 log 



A3 log 



(-) 



where the standard result for /i* is obtained, 

^J■c 



The coefficient oq is given in the appendix. This yields 
the non-trivial solutions 



Ai = 2Zujpiie 



' A-f,J(l + aoA) 



2ZwphAi* 



A 



(82) 



i + Mciog(:g^)- 



(83) 



Note that the gap in the spectral function Agp, found 
from the pole of the Green's function, occurs at j^l = 
Ai/Z in this approximation. 

Accounting for the approximations made in this deriva- 
tion we can use this form with three fitting parameters 
as in Eq. ([?7|.^^ 



(84) 



Asp = ClWphe ^-Mj(l + 03A) , 



where Ci > 0, C2 > 1. We can determine the parameters 
Ci , C2 by fitting to the numerical solution of Eq. ([71)) for 
/ic = in the regime < A < 0.5. Note that in contrast 
to the analytical solution the result of the numerical so- 
lution depends in general on the bandwidth W. It only 
becomes independent in the large bandwidth limiti^ 

For = 0, D = 8/cjph — 80, we find with ci = 1.7, 
C2 = 1.07 good agreement of the formula ([M| with the 
numerical results for A < 0.5 as can be seen in Fig. [T^l 
Note that for the Holstein model and larger values of 



A, the result in Eq. 



underestimates Agp as already 



pointed out by Allen and Dynes 
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FIG. 12: The spectral gap Asp as calculated from the numer- 
ical solution of Eq. (|74p as a function of A for D/uip\-i = 80 in 
comparison with the analytical form in Eq. (|84|1 with the the 
parameters ci = 1.7, ci = 1.07. 



(-)' 

(81) 



For fixed A = 0.5, we also compare the /ic-dependence of 
the analytical result in Eq. ([84p with the full numerical 
solution. In Fig. [T3| a comparison can be found, where 
we used the same values for ci and C2 as for /ic = and 
found that C3 = 0.8 gives a reasonable fit. 
We analyze the situation including only the C/^-term 
now. We set A3 = in this case and omit the con- 
stant f/-term in Eq. ([75)) . To determine the parame- 
ters Ai and A2, we use the following two conditions: 
S2i(0) = Ai = Y.^Nu/^,, and ^2i{iD) + S2i(-zD) ~ 
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FIG. 13: (Color online) The spectral gap Asp as calculated 
from the numerical solution of Eq. (|74p with different kernels 
as a function of for D/upi^ — 80 in comparison with the 
analytical result in Eq. (I84p . 



2A2/(1 + biD + b2D^) = J2t N'2^A,. This yields the self- 
consistency equation, 



l = A^ii 



N12N: 



21 



(85) 



Making similar approximation as in the first order calcu- 
lation we find for the coefficients Nij , 



and 



ith 



N2i = - 



1 



(A-aA*^)log(Ti), 



aiJ-c 



login 



Ti 

One also finds 



Zl + biD + b2D^ 
Ai 



(86) 



(87) 



N12 = ai2X - afil log f — A12 



and 



22 



l + biD + b2D^ 



log (—A22 



(89) 



(90) 



The expressions for ai2, A12 and for A22 can be found 
in the appendix. The coefficients A12 < 1 and ^22 < 1 
appear due to additional terms to 1/uj in the integrand 
which accelerate the decay towards higher energy. Hence, 
they lead to a reduced effective bandwidth, D — >■ ^22^* 
as discussed iir the previous section. The solution of Eq. 
yields 



Ai = 2Zujpiie 



)1) 



A 



22 



(92) 



and 



1 



log 1:^^22 



(93) 



The result for /i* has the same form as what was derived 
in Eq. (|73p . In addition a higher order term fi* i ap- 
pears. In the first order calculatioirs all higher order con- 
tributions cancel in the numerator of the expression for 
/i*. However, in the second order calculation this is not 
the case anymore and an additional term remains. The 
coefficient log(^22/^i2) does not increase with the band- 
width W and therefore the whole term becomes small in 
the large bandwidth limit due to the logarithm in the de- 
nominator. For the relevant values one has ^22/^12 ~ 2 
such that the coefficient of the /i^ is relatively small and 
the term does not contribute much for small values of ^c- 
For larger values, however, it does play a role. Thus it 
can account for the discrepancy between the irumerical 
and the analytical result from Eq. ([75)1 observed in Fig. 
[TUl where the numerical result does not seem to saturate. 

We would like to check these analytical findings with 
the full numerical solution of the self-consistency equa- 
tion. Formally, we had found very similar results to 
the ones of the last section derived from the projection 
scheme. For the comparison with the numerics we take 
these expressions and use the values for the parameters 
derived there, i.e. we use ^22 — ^ ^2 , where A2 = 0.197 
or log(/l22) = —1.624. We omit the term n* i, and for ci, 
C2 we take the same values as above, and we use 012 = C3. 
The reason for this procedure is that due to a number of 
approximations involved in the analytic calculation the 
results for these coefficient do not tend to be very accu- 
rate. Moreover, we aim for a unified description with as 
little parameters as possible. 

The result for the "only J7^"-calculation is added in 
Fig. [T3] and compared with the full numerical solution. 
We find good agreement. For small values of /ic the re- 
duction of Agp is smaller than for the first order term, 
but then Agp drops more rapidly. Notice that the second 
order result for fi* is analogous in the form to the first 
order calculation. The difference is the factor A22 < 1 in 
the logarithm in the denominator. Hence, the retarda- 
tion effects are less effective in this case as discussed in 
the previous section. The results here are derived inde- 
pendently from the arguments of the last section, but are 
fully consistent with them. The higher order term ^* is 
not very important for the values of /ic appearing iir Fig. 

m 

In the general case iircluding both the first aird second 
order terms in Eq. (j75p . we have to solve for three pa- 
rameters and hence have three self-consistency equations 
to be solved. In general this can be written as a matrix 
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equation A = MA. Algebraically this becomes rather 
lengthy and yields a number of different terms for n* , as 
discussed in the appendix. To simplify the discussion, we 
use results in the form of Eq. ([5^ with /i* as introduced 
in Eq. (EH), 



Mclog(^) + aM?log(^A22) 



1 



(94) 



for comparison with the numerical results. In Fig. I13l we 
have included the numerical result of Eq. (|74l) with the 
full kernel in Eq. ([75|) . We also included the analytical 
description based on Eq. ([M)) and ([7^ with the same 
value log(A22) = -2.28 or A22 = 0.1022 as in the pre- 
vious section. We find a rather good agreement for the 
range of values of fj,c- 

We can accurately calculate the dependence of the 
spectral gap Agp on fic for given A, D and Wph numeri- 
cally. However, it is not possible to calculate /i* directly 
from the self-consistency equation (|74|) . If we assume 
that the form of Eq. (|M)) . which neglects higher order 
terms of the form Eq. gives a good description and 
we can solve for /i*, 

ZC2 



A 

I + C3A 



log 



)(1 + C3A)' 



(95) 



As higher order terms are neglected this is not the com- 
plete result for fi* for the whole range of fic. However, if 
fj,* is interpreted as the quantity in competition with A 
to cause superconductivity, then this form is useful and 
the inversion of Eq. ([M)) can give us an estimate for /i* . 

The results for /i* obtained from Eq. ([M]) together 
with the analytical estimates in Eqs. (|ll73l72p can be 
found in Fig. [HI 
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FIG. 14: (Color online) /i* as a function of for D/ujph = 80 
computed from Eq. H95|l . We show in comparison the result 
for the first order and the second order calculation. The fit 
parameters are ci, C2 and C3 are the same as above. 

The results look very similar to the ones in the previ- 
ous section in Fig. [S] For small /ic, fi* increases linearly 



with fic in Eqs. (Ill72p . This implies an initially rela- 
tively large drop of superconductivity once the Coulomb 
repulsion becomes finite (see Fig. [T5)) . However, then the 
curves bend and the analytical results saturate. Only for 
very small values of fic the first and second order results 
are in agreement, otherwise the second order result is 
substantially larger. For the first order calculation the 
upper boundary is given by /^c,max — ^/^'^SiTJY^ inde- 
pendent of fic- It is reached both in the limit of large fic 
and large ratio Z^/wph. For instance for D/ujp\^ — 1000, 
we have fi* « 0.145, similar to what is usually used in 
the literature. Using fic — 0.5 for the screened Coulomb 
interaction gives with the same ratio for D/wph the value 
fi*c =0.112 (see Tab. U). 
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TABLE I: Exemplary values for fi^ for fic — 0.5. 
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ft* in Eq. ((44| 


0.366 


0.199 


0.136 


fi*c in Eq. ((72j 


0.700 


0.268 


0.166 



TABLE II: Exemplary values for fi* for fic 



We can see that the higher order results for /i* are 
generally substantially larger than in the first order cal- 
culation. In particular the results are larger than the 
simple minded estimate, ft* in Eq. (j44p . where retarda- 
tion effects are the same for the first and second order 
term. For this expression the limits of large bandwidth 
and large fic lead to the same result /I* = 1/ log(;j^), 
which is however reached already for smaller values, e.g. 
for fic = 0.5 we have /l* = 0.1233. In contrast, the re- 
suh in Eq. ^ goes to = l/log(^A22) in the 

limit of large fic and to fi* ^^^^ = 1 / \og{-^) in the limit 
of a very large bandwidth. With the estimate for A22 
above we find for D/ujp\^ = 1000 that in the large fic 
limit fi* jjjg^^ ss 0.216 is about 50 % larger than the first 
order estimate. For comparison we give some results in 
tables H] and [H] for fic — 0.5, fic = 1, and the ratios 
D/ojph = 10,100,1000. Notice, that the more accurate 
result for fi* as obtained from the numerical calculation 
and shown in Fig. IHlcan be still significantly larger than 
the estimate in Eq. ([7^ . We conclude that the usual 
result in Eq. ([T]) substantially underestimates fi* for in- 
termediate and larger values of fic- However, for very 
large values oi D / ojph retardation effects are operative in 
all cases and lead to a strongly reduced value of /^* . 
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C. Calculations for Ta 

For completeness we also include a brief section on the 
critical temperature Tj,. It is analyzed similar to the last 
section. The basis for the calculations is the pairing ma- 
trix in Eq. . For the pairing vertex we use the same 
terms as in the previous section in Eq. (|75l) . 



(96) 



The effect of the Z-factor is neglected. Then with the 
bare Green's function and a semi-elliptic DOS 

X (^cJ„l;0,0) = = pn— Wl + — - 1 . 

(97) 

n(zw„) is calculated numerically from the free Green's 
function. With this we compute the matrix M^-^ n^ in Eq. 

and search for the largest eigenvalue. We compare 
the results for Tc for three different calculations: (1) in- 
cluding only the [/-term, (2) including only the C/^-term 
and (3) including both. 

Apart from the first order calculation it is not easy 
to find a good analytic approximation for the eigenvalue 
equation p2p . In principle, one can do something similar 
to what has been done in the last section and make an 
appropriate ansatz for the eigenvector. To first order in 
U this works reasonably well. We find a result of the 
standard form, 



Tc = ciWphC ^-f;<i+-3A) 



(98) 



where fi* is given by Eq. ([85)1 . For the higher order 
analysis we did not pursue an analytical solution, and 
instead also assume Tc as in Eq. ([55]) and for fi* the 
form in Eq. (|72|). 

First we fix the constants ci and C2 by fitting to the 



result for fic — 0, see Fie 



We take the value ci = 



1/1.2 « 0.833 as in Ref. |5J and find a good fit for C2 = 
1.04. Similar as before the agreement is only good up to 
values of A < 0.5. 

In Fig. [in] we give the numerical result for Tc as a func- 
tion of (Uc for the first order calculation, only the second 
order and first plus second order. We have kept A = 0.5 
constant. 

Tc decays in a very similar way as the spectral gap when 
He is increased. We included results from the analytic 
fit form in Eq. ([M[) with the value C3 = 0.8. For the 
first order calculation, we find good agreement with the 
numerical result for the case A = 0.5. For the second 
order results we use the same parameters as before and 
the fits are reasonable, i.e. they show that the functional 
form is close to the actual result. 

Note, that in spite of the identical form for the pairing 
kernel, the results for Tc here and the ones for the gap at 
T = in the previous section are obtained from two in- 
dependent calculations. The results indicate that for the 
range of parameters studied the analytical expressions 
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FIG. 15: (Color online) Tc as a function of A for D/wph ~ 80. 
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FIG. 16: (Color online) Tc as a function of /ic for D/ixiph — 



Eqs. (|ll73l72p describe the effect of /i* on Tc quite accu- 
rately. In particular the results are very similar to what 
has been found in the projective approach before, and 
hence a consistent picture emerges. We conclude that 
higher order dynamic effects, and in particular the sec- 
ond order contributions, give an important correction to 
the usual results for the Coulomb pseudopotential. Not 
only does the higher order term give a direct increase of 
the coupling, as would be the case in expression in Eq. 
([44| it also leads to a reduced effective bandwidth in the 
logarithm in the denominator. This is an effect which to 
our knowledge has not been discussed in the literature so 
far. In the following section we analyze how these effects 
are manifest in the non-perturbative DMFT calculations. 



VI. DMFT AND PERTURBATIVE RESULTS 

In this section we put the analysis of the previous 
sections together and compare with non-perturbative 
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DMFT calculations, which include all possible renormal- 
ization effects. We first clarify that it is therefore very 
important to work with renormalized parameters for the 
interpretation of the results. Then we demonstrate up 
to which interaction strengths the perturbative results 
derived in the previous section are reliable. 

We focus on calculations for the spectral gap Agp at 
T — and consider a half filled band. The spectral gap 
Asp is extracted directly from the gap edges of the di- 
agonal spectral function. It is usually well approximated 
by the product zE2i(0), but it can be a bit larger due to 
the frequency dependence of the off-diagonal self-energy. 
For an interpretation of the DMFT results we need to 
compare with the perturbation theory (PT) results. We 
include the following terms for the self-consistent PT cal- 
culation: For the diagonal and off-diagonal self-energy we 
use Eqs. ([5T]) and ([55)) . The vertex is approximated by 
the contributions from RPA screening in Eq. ([M)) up to 
second order in U and the second order in U corrections 
in Fig. [S] For a small phonon frequency, the uj — value 
of the vertex function provides a good approximation. 
From the diagrams in U we take Eqs. (|^^ . and ([M]) - 
(j57p into account. The phonon propagator is taken as an 
input from DMFT calculations and not calculated self- 
consistently, similar to what has been done in Ref. [ij. 

To get a feeling for the renormalization effects let us 
first consider calculations, where the bare parameters 
•^0 = Po2(7^/wo and ujq in the Hamiltonian Q are kept 
fixed. For /ic = the system has a superconducting solu- 
tion and we analyze how the gap Agp is affected, when fic 
becomes finite. For Aq = 0.308 and Aq = 0.382, the cor- 
responding spectral gaps Asp are shown in Fig. [T71 For 
comparison we have also included the results from the 
self-consistent PT as explained above in Fig. [T71 These 
are seen to be in very good agreement with the DMFT 
result. 




FIG. 17: (Color online) Behavior of the spectral gap Asp and 
zE°"(0) for constant Ao = 0.308 and Ao = 0.382 {ujo = 0.1 in 
both cases) as a function of fic- 

We find a rapid decrease of Agp when fj,c becomes finite. 



similar to the results in Fig. [T31 However, notice that 
Asp goes to zero for even smaller values of fj,c ~ 0.3. 
Naively one could conclude that fj,* is even larger than 
what was discussed in the last section even though we 
are still at very weak coupling in [/. As we will see, how- 
ever, this drastic reduction of Asp has a different origin. 
It can be understood by analyzing the relevant effective 
parameters and the PT results. 

Since, as shown in Sec. IVl manv aspects of the PT are 
well described by approximate analytical results, such as 
Eq. (|7^ for ^* or Eq. for Agp, it makes sense to 

use those equations to analyze the results. In Sec. IVlthe 
equations were used in terms of bare parameters, but to 
understand the DMFT results we need to use renormal- 
ized parameters. Those effective parameters extracted 
from the DMFT calculation are shown in Fig. [18] as func- 
tion of flc- 
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FIG. 18: (Color online) The effective parameters (a) the 
quasiparticle renormalization z, (b) the renormalized phonon 
frequency uq/loq (c) the renormalized coupling constant 
{g^ I gf' (the full line corresponds to the result with a and h 
obtained from the free Green's functions) and (c) the effective 
A as a function of \ic- 

The 2:-factor in Fig.[TH](a) increases moderately with /ic- 
According to Eq. ([M[l with Z = this would actually 
help superconductivity, so it can not be responsible for 
the reduction in Fig. 1171 The renormalized phonon fre- 
quency ujq in Fig. [15] (b) also increases with fic showing 
that corrections to the phonon self-energy and electron- 
phonon vertex from finite U are important. In the pref- 
actor of Eq. (I84p with Wph — (j-!q this leads to an en- 
hancement of Agp, whereas fj,* increases with loq which 
leads to a reduction of Asp. Both effects are not strong 
enough to be decisive. The strongest and most impor- 
tant effect of flc can be seen in Fig. [TS] (c) and (d) , where 
we plot the renormalized coupling from Eq. (j39p and 
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the effective A according to Eq. The parameters 

a and b in Eq. ((5^ are calculated from the full Green's 
functions. We find that the renormalized coupling de- 
creases substantially with /ic. In addition the phonon 
self-energy is modified for finite U. Therefore, the effec- 
tive A becomes much smaller. So the main reason for the 
rapid suppression of Agp is the strong effect of the cor- 
rections to the electron-phonon vertex and to the phonon 
self-energy, such that the effective A decreases. 

Using values Wph — uJq we can calculate results for the 
Coulomb pseudopotential /i* according to Eq. (|72|) . The 
ratio of electron over phonon scale is with D/ujq ^ 35 not 
as large as in the last section, and thus retardation effects 
not as effective. For the largest value fic — 0.318 we find 
/i* ~ 0.2. At this /ic we have A ~ 0.45 and the analytic 
expression in Eq. ([5^ yields Agp/wo ~ 0.0002, where 
Z = z^^, i-^ph = wj, and for the fitting parameters Ci 
take the same values as in Sec. IV Bl This is in agreement 
with the DMFT finding that superconductivity goes to 
zero then. If we compare the results for the spectral gap 
according to Eq. (|84l) with the renormalized parameters 
in Fig. [IS] with the results from the full calculation in 
Fig. [T7|we find good agreement. Hence we conclude that 
the superconducting state at small /ic can be well under- 
stood in terms of the effective parameters A, wj, z and 
fi* and the approximate equations for /i* and Agp. 

The conclusion up to this stage is that the strongest ef- 
fect of the Coulomb repulsion in the DMFT calculations 
is to renormalize the effective A via electron-phonon ver- 
tex and phonon propagator such that superconductivity 
drops to zero rapidly. Even at very weak coupling these 
effects play an important role in the Hubbard-Holstein 
model and must be taken into account. Since in this 
work our main interest is the effectiveness of retardation 
effects visible in the direct competition between A and 
/i*, we will offset the vertex-correction effect in the fol- 
lowing. We do this by appropriately adjusting the bare 
parameters, i.e., increasing Aq and cjq- We proceed by 
keeping A — A(/ic) as defined in Eq. ([^^ constant. For 

we rely on the perturbative results. This was found 
to give a relatively good description up to [/ ~ 
First we do a calculation for certain bare values ujq and 
Ao at /ic = 0. From the phonon spectral function we can 
then extract the value and A from (j29p where, g'' = g. 
Then we choose a finite value of /Xc- We have to in- 
crease Ao and wq, such that roughly equals the [7 = 
value and A remains approximately the same according 
to Eq. (P^)) . For the renormalized coupling g^ entering 
Eq. (|29p we consider two conditions: (a) the second or- 
der expansion as in Eq. (|4ip and (b) the RPA series plus 
the second order terms as in Eq. (l39l) . With the set of 
bare parameters (Ao,woiMc) found from this procedure 
we also do perturbative calculations for comparison. We 
distinguish PT-a, which includes up to second order di- 
agrams for the electron-phonon vertex, and PT-b which 
includes the RPA series plus the second order terms for 
the electron-phonon vertex. The results of such a cal- 
culation, where A ~ 1 according to condition (a) and 



loq/D ~ 0.025 are shown in Fig. [12] as Agp vs fic- 
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FIG. 19: (Color online) Behavior of the DMFT result for the 
spectral gap Agp and zT,°^{0) for constant A ~ 1 according to 
condition (a) for the vertex corrections (see text), ujq ~ 0.05 
as a function of /ic in comparison with PT-a and PT-b (see 
text) and the analytic formula in Eq. (|84|) . 

The DMFT results show a steady decrease of Agp on 
increasing /ic. Since the effective A and cjph are kept 
constant the diminution of Agp is now due to the com- 
petition with the Coulomb repulsion fic- To understand 
the result quantitatively we compare it with the pertur- 
bative calculations PT-a, PT-b and the analytic results 
based on Eq. ([M)) using the effective parameters and fi* 
according to Eq. ([7^ with the value of A22 as in Sec. |V] 
We find a relatively good agreement of the DMFT result 
with PT-a and the analytic formula up to fic ^ 0.5. This 
demonstrates that (i) the effective parameter description 
is appropriate, (ii) the electron-phonon vertex correction 
according to condition (a) is suitable and (iii) that the 
derived higher order form for the Coulomb pseudopoten- 
tial in Eq. f72|) captures correctly the results of the PT 
and the of the full DMFT calculation. This validates the 
analysis of the Sec.|V]in a more complete calculation and 
it corroborates our findings for /i* and retardation effects 
by comparison with non-perturbative DMFT up to inter- 
mediate values of fic ~ 0.5. For larger values of fic, higher 
order correction enhance the value of /x* such that Agp 
is suppressed stronger. In addition the electron-phonon 
vertex is not well described by condition (a) anymore. 
The PT-b calculation agrees with DMFT quite well up 
to values fic ^ 0.3, but then overestimates the reduction 
of the electron-phonon vertex and therefore leads to a too 
strong suppression of superconductivity. 

To extend the present analysis to larger values of fic, 
one needs to find reliable estimates for the full renor- 
malized electron-phonon vertex F^J''''. A higher order 
perturbative analysis or non-perturbative calculation for 
this quantity could provide this. A consistent calculation 
up to certain order in fic should then also include higher 
order corrections to the self-energies, which will compli- 
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cate the analytical calculation for /i* further. This is 
beyond the scope of this work. 



VII. CONCLUSIONS 

For the occurrence of conventional superconductivity 
it is important that on the one hand there is a sizeable 
electron-phonon coupling and on the other that the detri- 
mental effects of the Coulomb repulsion are reduced suf- 
ficiently by screening and retardation effects. We have 
stressed in this article that while the effect of the for- 
mer can be described in a controlled fashion by Migdal- 
Eliashberg theory, the standard approach to the latter is 
based on an uncontrolled approximation, since there is 
no Migdal theorem for the Coulomb interaction. For the 
Hubbard-Holstein model we have analyzed this issue in a 
controlled framework by a combination of perturbative 
calculations and non-perturbative DMFT calculations. 
We have shown that the conventional arguments based on 
the lowest order diagrams become modified when higher 
order corrections are taken into account. There is still 
a reduction of the Coulomb repulsion due to logarithmic 
terms in the denominator. This demonstrates that indeed 
a small value of fj,* can be obtained due to retardation ef- 
fects even when higher order corrections are considered. 
Thus our results support the arguments by Morel and 
Anderson qualitatively. However, the effective energy 
scale separation is reduced. This is due to the dynamic 
behavior of the higher order diagrams which make retar- 
dation effects less operative. This result was shown ex- 
plicitely in two independent calculations in a combination 
of analytical and numerical arguments. By studying the 
occurrence of superconductivity with DMFT, where all 
higher order corrections are present, we were able corrob- 
orate our findings up to intermediate coupling strength. 
The perturbative approach allowed us to distinguish dif- 
ferent renormalization effects. Our analysis is limited to 
intermediate coupling strength due to the difficulty to 
reliably estimate the vertex corrections to the electron- 
phonon vertex. Non-perturbative calculations for this 
quantity would be desirable. 

We conclude that the usual expression for /i* in Eq. ([T]) 
is only valid for small values for fic ■ The second order cor- 
rections lead to less efficient retardation effects as shown 
inEq. (1721), 

= r^^^ 7 (99) 

where we found a ~ 0.1 for a typical large energy separa- 
tion of electron and phonon scale i?oi/wph. The coefhcient 
a is given by the limit — > of the particle- hole bubble 11 
divided by po, the DOS at the Fermi energy, and a can be 
estimated from the decay of H with uj. In addition higher 
order terms appear such that /i* does not saturate in the 
limit of large fic- Thus, for systems with sizeable effective 



Coulomb repulsion we should expect larger values for fi* 
than what is traditionally quoted. 

It is premature to draw detailed conclusion from our 
calculations for real materials such as lithium, since the 
Hubbard model is not an accurate description for itin- 
erant metallic systems, which are better described by an 
electron gas model. Nevertheless one can understand our 
results as a qualitative trend, which shows that for sys- 
tem with less efficient screening, such that is larger, 
we expect an enhanced value of p,* as compared to the 
traditional quote /i* ^ 0.1. This could be expected in 
systems with larger values of r^. Hence our conclusions 
would seem to be line with the observation in a number 
of of materials, such as Li with enhanced /i*, and what 
is quoted in work based on DFT calculations. Similar 
conclusions are important for systems where the ratio 
Eci/idph is reduced from the usual scenario, which could 
be the case in picene. We would like to stress, however, 
that a more accurate attempt of understanding the prob- 
lem should also include the dynamic effect of screening 
the bare Coulomb repulsion in a metal, which was not 
taken into account explicitely. As discussed in the in- 
troduction this can lead to very small and even negative 
values of /u*. We expect that the combination of these 
effects and the corrections studied here, which lead to 
an enhancement, will eventually lead to the physical val- 
ues occurring in nature. More detailed calculations are 
required to fully resolve this quantitatively. 
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Appendix A: Details for the calculation of the 
spectral gap 

In this appendix we collect some details for the analytic 
calculation in Sec. IV Bl 

1. Integrals for the first order case 

For T,2i{0) we approximate the integrals over small fre- 
quencies, 

^ J duj G2i{iuj)g^D{-iuj) ~ ^ log(ri) , (Al) 
and 

^ i do. G2i{iuj)U = ^ log(ri) , (A2) 

-"ph 
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where Ti is given in (j88p . For larger frequencies we write 



1 

2^ 

and 
1 

2^ 



h D 

duj + J dw 

D liJph 
-"ph D 

duj + J du! 



G2i{iuj)g^D{^iio) = aoAsX , (A3) 



G2i{iuj)U = -Aa^clog 



We have 



D 



Oo = / da; ^ = - log | 2 



W 1 + ii^- 2 

ph Pli 



D2 



(A4) 



(A5) 



For E2i(i£'), the electron-phonon contribution is small 
for ujph_/D <C 1, and the Coulomb contribution is the 
same as above. This yields Eq. (|5T|) . 



The integral on the left hand side can be carried out an- 
alytically, but the expression is lengthy and not instruc- 
tive. We can express A22 as 



A22 — exp 



l( h{D) / 1 

u\2h{u)\h{D 



+ w) hiD-u) 



- 1 



Since the function l//i(_D + cj) + 1/ fi{D — uj) increases 
in the integration interval the coefficient A22 comes out 
larger than A12 . For the parameters above we find A22 ~ 
0.44. The coefficient for A reads 



0-12 



did 



UJ 1 



(A9) 



Wph 



The integral can be solved analytically but the expression 
is lengthy. Due to the reduction factor l//i(a;), ai2 term 
is a bit smaller than qq in Eq. (|A5p . 



2. Integrals for the only second order case 

We give some results for the second order calculation. 
To determine the coefficient A12 we calculate 



where 



fl{uj) = 1 + bi\u\ + b2UJ^ . 



(A7) 



One finds that 



A, 



exp 



1 (2b2-{hl + bib2D) 



V MD) 

2^2- (bf + bib2^ph)- 
/i('^ph) 



) exp(^ 



h(.D) 
bl - 66162 



fa 



with i?i = -^Z 462 — 61 and 
'61+262!) 



/a = arctan 



-arctan 



61 + 262a;ph 



(A8) 



By comparing the integrals in Eq. (IA4[) and Eq. (|A6P 
we see that the factor < ^12 < 1 leads to a reduction 
of the effective bandwidth. For the given fitting values 
61, 62, and large ratios D/uph > 100 one finds A12 ~ 0.2, 
which can be compared to the calculation in Eq. ([7T|) . 
The coefficient A22 is obtained from the integral 



D 



1 



+ uj) fi{D-uj). 



2l0g(^A22) 



3. Calculation up to second order 

We use the following three conditions: (i) S2i(0) = Ai, 
(ii), 

S2i(i6i) + S2i(-i6i) ~ 2A3 -t- 2A2//i(6i), (AlO) 
with 61 = 1/61, and (iii), 

J:2l{^D) + I]2i(-*i?) ^ 2A3 + 2^2/ h{D). (All) 
We use I]2i(0) =E»^hA^, 

S2l(»6i) + S2l(-»6i) = Y.^ N2^\, (A12) 

E2l(^i?) + S2l(-^i^) -Ez^3»A,. (A13) 

This implies Nu — AIu, 

M21 = /i(6i)iV2i/2, M22 = ./i(6i)A^22/2, (A14) 
M23 = /i(6i)(iV23 - 2)/2, M31 = iV3i/2, (A15) 
M32 = (iV32 - 2/h{D))/2, Af33 = N33/2. (A16) 

The calculations give with certain approximations in the 
integrals 

Nn = |(A - Mc - a^c) log(Ti), (A17) 



N12 = ai2 A - He log ( A\y ) - afil log (— ^^2^ 1 , 



(A18) 



A^i3 = -A - log ( — ) - a^^ log (^^la^) > (A19) 



(- 



22 



N,, = - 2m.. log (—4?) - log ( — 4f , . 

(A21) 



D 



032 ^ I duj — 



, (A29) 



^23 = a.3A - 2^. log (^) - 1^ log (^4^3^ > 

(A22) 



033 



duj 



ph 



-), (A30) 



N3i = - 



A^32 = a32A - 2^c log ( A. 



(A24) 



1(2; 



^22 = GXP 



1 



1 



iJ\2h{u) \fi{bl+Lj) fi{bi-uj) 

(A31) 



iV33 = a33A - 2^. log ( — ) - 1^ log (^4f 
^Wph/ fi[D) Vwph > 

(A25) 

The coefficients 012, ^12 — ^12 and = A22 were 
given above. The others read, 



^23^ = exp 



1 



-1 



+ 0;) 

(A32) 



^(1) _ ^(1) _ 4(1) _ ^(2) _ ^ / AKh) 

^12 —^32 —^22 —^13 ^ y 1 



a22 = / do; 



1 1 



a23 = / dw 



Si 



/a 



(A26) 



1 1 

+ 



^pn ph ph 



(A27) 



(A28) 



A 



exp 



dcif/lMf^ 
c^V 2 V/ip 



1 



-1 



(A33) 

One can give analytic expressions for the integrals, which 
are however, lengthy and not very instructive. For typical 
values for ZJ/wph, we have A^^^ w 1.6 and ^gg-* sa 1.7. All 
other coefficients obey < A^""^ < 1, leading to a reduce 
effective bandwidth as discussed before. The coefficients 
small for ZJ/cjph ^ 1. 

From the matrix equation A = MA, we can derive 
the self-consistency equation. 



1 = Mil + + I M12 + - 

I-M33 V 1 



M32 \ M2i(l-A/33)+M23M3i 



- A/33/ (1 - M22)(l - M33) - M23M32 



or when using the expression for Mn with Ti = ^"p''+V(^|^ph) +4 i as in Eq. 



log(ri) 



1-M: 



33 



M12 



M: 



32 



Af2l(l-Af33) + M23M3i 



1 - A/33 / log(ri) (1 - A/22)(l - A/33) - A/23A/32 
I 



(A34) 



(A35) 



When solved for gap parameter Ai, the term in the square brackets is the exponent, such that terms to the 
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right of A contribute to the expression for the pseudopo- 
tential /i* . We now would hke to argue that in an expan- 
sion in /ic the dominant term is of the form of Eq. ([72]) . 
We neglect the terms involving A in Mij to simplify the 
arguments. We find then that — /ic — together with 
13 3i_iog(T]j_ gj^^gg ^ tei-jn of the form in Eq. ((7^ . 



, lO! 



V^ph 



^33 



(A36) 



where 



T(2) 
^33 



exp 



-1 



(A37) 



There are also additional terms to order /ij? and /i^ whose 
coefficients are not proportional to log(D/a;ph) in the nu- 
merator (cf. discussion in Sec. lVBl) . The denominator in 

the other part, F2 = {I - Af22)(l - M33) - M23M32, 



has similar properties to the one in Eq. (jA36l) but con- 
tains also contributions to order /i^ and /i^ and terms 
^ log(L)/ajph)^. Due to a cancellation the lowest or- 
der term in the numerator, M2i(l — M^-^) + M2iM^i, 
is ^ a/i^(/i(6i)//i(Z?) — 1). From the prefactor, 



Afi2(l-M33)-HM; 



32 



l-M: 



33 



(ASS) 



the lowest order term is — l/[/i(Z3)(l — Af33)]. This gives 
a contribution 



\fiiD) 



1 



h{D)F2 



(A39) 



which is smaller compared to a/x^. All other terms are of 
the order jjL^ and higher. Hence, in an expansion in /ic 
the term to the right in the square brackets in Eq. (jA36P 
gives a smaller contribution to /i* . This explains why the 
numerical results in Sec. IV Bl were fit well by an expres- 
sion involving /i* of the form in Eq. (|72l) . 
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